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Abstract — A new design method for high rate, fully diverse 
('spherical') space frequency codes for MIMO-OFDM systems 
is proposed, which works for arbitrary numbers of antennas 
and subcarriers. The construction exploits a differential geomet- 
ric connection between spherical codes and space time codes. 
The former are well studied e.g. in the context of optimal 
sequence design in CDMA systems, while the latter serve as 
basic building blocks for space frequency codes. In addition a 
decoding algorithm with moderate complexity is presented. This 
is achieved by a lattice based construction of spherical codes, 
which permits lattice decoding algorithms and thus offers a 
substantial reduction of complexity. 



I. Introduction 

In MIMO (Multiple Input Multiple Output) systems space 
time coding schemes have been proven to be an appropriate 
tool to combat fading holes of Rayleigh flat fading channels 
and exploit the spatial diversity gains. For frequency selective 
SISO (Single Input Single Output) channels (due to multipath 
propagation) OFDM (Orthogonal Frequency Division Multi- 
plexing) provides an efficient strategy to obtain a collection of 
mutually orthogonal flat fading channels on each subcarrier. 
The number K of subcarriers in such a system is not a free 
design parameter but predetermined by the bandwidth (which 
determines the sampling rate) and the coherence length of the 
channel (which determines the maximal OFDM symbol length, 
which is the maximum number of subcarriers). MIMO-OFDM 
is the combination of MIMO techniques to frequency selective 
Rayleigh (multipath) channels. Space frequency coding aims 
to combine the advantages of space time coding with the 
multipath diversity from each channel tap. 

A corresponding code design for fully diverse space fre- 
quency codes C has been proposed in [1]. In their work 
they separated the code construction into the design of an 
appropriate transition matrix A incorporating the multipath 
diversity and the construction of a corresponding space time 
code C. This method has the advantage of splitting up the 
design problem into two separated problems by keeping the 
multipath echoes of the transmit signal mutually orthogonal. 
For these reasons this work focuses on the design method [1], 
although other approaches have been investigated (e.g. [2]). 

To obtain maximum rates by retaining full (spatial and 
frequency) diversity order the block length T of the 'inner' 
space time code has to be as large as possible (i.e. T — K/L, 
whereas L denotes the number of taps). This requirement rules 
out most of the space time code designs considered in the 
literature so far. Common research has been focused mainly on 
small dimensional coding spaces (i.e. T = n t ) whose elements 
consist of square n t x n t matrices, whereas nt denotes the 



number of transmit antennas. These designs correspond to the 
case K/L = n t only. To attack the more general (space time) 
code design problem an inspection of the maximum likelihood 
decision criterion reveals, that it is a good first approach to 
consider code design as a packing problem. The corresponding 
coding space turns out to be the complex Stiefel manifold of 
complex T x n t matrices with unit norm mutually complex 
orthogonal column vectors. 

On the other hand the design of optimal (Welch bound 
equality WBE) signature sequences in CDMA (Code Division 
Multiple Access) leads naturally to a packing problem on 
spheres [3], [4]. Explicit constructions of spherical packings 
can be found on N.J. A. Sloane's web page 1 . 

The first achievement of this work is to employ differential 
geometric methods to establish a connection between spherical 
packings/codes and (general) space time packings/codes, such 
that the rich variety of already known spherical codes can 
be transformed to space time codes in various dimensions. 
This mapping preserves (geodesic) distances, and therefore the 
performance at the receiver can be controlled (up to a certain 
degree). The so constructed codes are called spherical space 
time codes and the resulting space frequency codes (obtained 
from the procedure described above) are named spherical 
space frequency codes from now on. 

The second part of this work is concerned with decoding 
strategies. The need for high rate code designs demands 
advanced decoding algorithms. Even moderate numbers K 
of subcarriers need huge codebooks of size 2 KR in order to 
retain a specified rate R, making brute force maximum likeli- 
hood decoding impossible for practical applications. Based on 
wrapped spherical codes [5] which are decodable in terms of 
lattice algorithms, the connection to space time codes will be 
exploited to do the same for space time and space frequency 
codes. This preferred decoding strategy prohibits high rate 
spherical code classes obtained elsewhere, such as Sloane's 
dense packings or packings obtained from frame theory (e.g. 
[6]). 

On the other hand wrapped spherical coding offers a wide 
range of lattice based constructions which can be adapted to 
arbitrary numbers of antennas and subcarriers. Moreover for 
each lattice of the appropriate dimension a target minimal 
distance can be prescribed, leading to a corresponding wrapped 
spherical code which can be processed further to obtain 
space time and space frequency codes. In this way whole 
(fully diverse) families of codes can be constructed for each 
prescribed number of antennas and subcarriers, whereas the 
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rate can be adapted by choosing appropriate minimal distances. 

The text is structured as follows. In the next two sections 
the channel model and the basic code design criteria will be 
derived. Section IV is devoted to the design of space frequency 
codes, which has been made popular in [1]. Section V sketches 
the differential geometric connection between spherical codes 
and space time codes and describes the corresponding map- 
ping. In section VI the second part of this work starts. It 
contains the encoding and decoding procedure (based on 
[5] for the spherical part) for the proposed spherical space 
frequency codes, which are simulated in section VII. Finally 
in the last section the results are summarized for a quick 
reference. 

II. Channel model 

Let us adopt the following notations: x = FFT(a;), x' = 
IFFT(:e) = FFT^(x), 7 denotes complex conjugation, •* the 
transpose, -t the hermitian transpose, 1 denotes the identity 
matrix and the all-zero matrix is denoted as 0, whereas their 
dimensions will be clear from the context, if not indicated 
by some subscripts. Using the common (complex conjugated) 
K-th root of unity u> = the (if -point) (I)FFT matrix 

reads (i)fft = [(l)FFT fefe /]£~;^ with components FFT fefe / = 
, resp. iFFTfefc/ = k . The frequently used 

operations of (frequency) modulation p, = diag(w /c )^Z 1 and 
(time) translation r = [e^+i) moiK] k =a on ^> K ( e fc denoting 
the fc-th standard basis vector) transform into each other by 

H = FFT°T°IFFT, resp. T = IFFT°/X°FFT. 

The channel model assumed here is the frequency selective 
Rayleigh fading MIMO-OFDM channel with K subcarriers, 
L taps Hi G C" tX ™- ~ i.i.d. Cff(0,±), and n t , (resp. 
n r ) transmit (resp. receive) antennas. Before transmission a 
K point IFFT is applied to each K x n t MIMO-OFDM 
code symbol C = (cq, . . . , c^-_ 1 ) t , which is then prepended 
by the cyclic prefix of length L. The receiver applies a K 
point FFT and removes the cyclic prefix before (maximum 
likelihood) decoding. Then, the transmission equation reads 
in the frequency domain 



Cfc 



yfpc k H k +n k , k = 0, . . . , K - 1 



whereas 



L-l 



(i) 



(2) 



and n k € C lxn *' (noise), with both H k and n k having i.i.d. 
CAf(0, 1) distributed components for all k = 0, . . . , K—l. It is 
assumed that the transmitter knows only the channel statistic 
while the receiver perfectly knows the channel coefficients. 
The rows c k of C are assumed to have unit norm, such that p 
is the SNR (signal to noise ratio) per subcarrier at each receive 
antenna. 

III. Space frequency code design criteria 

This section describes the code design criteria developed in 
[1] with emphasis on the underlying geometrical framework. 



Let us express the transmitted symbols C as C := sj p^ 
assuming $ to be an element of the complex Stiefel manifold 



Vn t , K = {* e € Kxnt I = 1} 



(3) 



which is the basic coding space from now on, compare [7],[1]. 
A space frequency code is given by a discrete set C C V® k 

and provided the received signal C — p-^ + N the 
maximum likelihood decision reads (see [8]) 



ML 



arg mm 




(4) 



In the sequel let us adopt the notation that in order to 
emphasize the space Frequency (code) or space Time (code) 
interpretation of a given quantity subscripts - F , respectively 
■ T will be added. Further, for the difference symbol A = 
$ — 5' of two codesymbols we define the 'multipath extension' 
A mp = (A, pA, p L ~ 1 A). Then the Chernov bound for the 
probability of mistaking $ for * is given by (compare [1]) 



/n t L 



ch F = n t 1 + <?f°?(aj] 



(5) 



whereas g F := jp-^j; and a(M) = (o-j(M)) generically 
denotes the vector of singular values of the matrix M in 
decreasing order. Expanding the product leads to the diversity, 
resp. diversity sum d F and product p F 

n t L n t L 

Vzv F := J] [1 + g F aUAJ] = ]T (6) 

i=l i=0 

whereas (sym, denotes the z-th elementary symmetric poly- 
nomial) Sl := sy mi (CT?(A mp ), . . . , al tL (A mp )) and 



K ■■= si 

2 



I A. 



:= s ntL = det(At p A mp ) < (p 



i||A|| F ^ = id^ 

2 \L 



(7) 
(8) 



with d T := || A || F , p T := V dct At A denoting the correspond- 
ing diversity sum and diversity product for space-time codes. 
Note, that due to d 2 = Ld 2 the maximum likelihood detection 

' F T 

for space frequency codes is essentially the same than for 
space time codes. Further, equality in (Fischer's inequality) 
(8) holds, if AV A = 0Vi=i,..., L -i, provided p T ^ 0. 
Thus a sufficient condition for the maximization of diversity 
(sum/product) is 



max : d 2 = 11$ - <J/|| 2 > 
max : p 2 = dct A* A > 



V i= i,...,L-i(At A1 ( A = 



A'tr'A' = 0) 



(9a) 

(9b) 
(9c) 



performed over all pairs $, * of code symbols satisfying 
$ ,\E' € Vl c r K . So, while the first two conditions are 
borrowed from the theory of space time codes, the third set of 
conditions is a serious constraint: The difference vector A of 
two code symbols $, $ has to be perpendicular to its row-wise 
l-th order modulated (resp. its Z-th order down-shifted) version 
for I = 1, . . . , L — 1, reflecting the strong correlation between 
each l-th channel realization in the frequency domain. The 



demanded orthogonality forces the best possible (frequency) 
diversity induced by the multipath delay spread. 

The homogeneity of the diversity (6) permits the interpreta- 
tion of (9) as a constraint packing problem on the Stiefel man- 
ifold V„ c 

l k- Fully diverse space frequency coding becomes 
a packing problem in V® t K (due to (9a)) subject to the 'space 
time' constraint (9b) and the 'space frequency' constraint (9c). 
For space time codes the corresponding interpretation has been 
pointed out in [9]. 

IV. Space frequency codes from space time codes 

As has been pointed out in [1], the design of fully diverse 
space frequency codes splits into two independent parts, 
reducing the overall complexity of the problem considerably. 
First construct a space time code C C V^.t (subject to the 
first two conditions in (9)) for some T. Then find some matrix 
A' satisfying 



A' =(A',TA',...,T L - 1 A>)eVV 



(10) 



The final space frequency code fulfilling (9c) is then given as 
[1] 

C = ¥¥T(A'C') (11) 

From this construction the rate R of the space frequency code 
(neglecting the cyclic prefix) and the rate R of the internal 
space time code are related by 



T - 
R=-R 



(12) 



from which it is obvious that high rate codes require large 
block lengths T. On the other hand to retain the maximal 
diversity order Ln t n r of the space frequency code T < 
K/L is demanded. Therefore the code design aims at code 
constructions satisfying 



(13) 



and we concentrate on this case in the sequel. 

As a further notational convention we introduce the tilde ~ 
for quantities corresponding to the 'inner' space time code C. 
Recall that C in (11) has full (frequency and spatial) diversity, 
if C is a fully diverse space time code and for large block 
lengths this is a complex code design problem. Therefore the 
strategy in this work is to concentrate on the simpler first 
order design criteria (9a) in the next section and develop some 
adjustments for full diversity later on (section VI), compare 
the remark in the last paragraph in section III. Due to the 
factorization (11) the goal investigated in the next section is 
to find dense packings in the metric space (V^ T ,dr)- 

Let us close this section by some analysis of the frequency 
diversity generating matrix A'. Random constructions have 
been proposed in [1], but here we present a deterministic 
construction which in addition is especially designed to benefit 
from a low crest-factor: Set 

A = I <-xi> | - ) ] (14) 



K 



k=0 



(it is well-known that this sequence belongs to the set of 
so-called bi-unimodolar sequences, i.e their FFT is as well 



unimodular). Then the matrix U(c) = (uo, . . . , Uk-i) '■= 
(c, tc, . . . , t k ^ 1 c) with c = FFT 1 ' (-^= A) is a circulant matrix 
which is also unitary. Let us define 



A' := {u Q ,u L , . . .,U( T -i)l) 



(15) 



then t 1 A' = (u u u L+u ■ . . ,U(t-i)z,+z) and A'^A'^ = 1 
as desired. This choice of A distributes the transmit energy 
equally in the time domain (compare [10] and references 
therein for an account to the construction of sequences with 
zero autocorrelation sidelobes and zero crosscorrelation). In 
general this construction scheme provides a real K dimen- 
sional parameter space A^ = e l<t>k , < (f>k < 2tt, k = 
0, . . . , K - 1 for A'. 

V. Space time packings from spherical codes 



The Stiefel manifold is canonically a homogeneous space 
(U(K) -^U V® >K ; U(K-nt)): The canonical left multiplica- 
tion of rij-frames in <C K by unitary KxK matrices transforms 
each pair of n t -frames isometrically (with respect to the metric 
induced from the unitary group) into each other. Thus the 
action of the unitary group U(K) on V® K is transitive and 
establishes the canonical diffeomorphism 



(16) 



Homogeneity is the natural generalization of spherical sym- 
metry: At any point the manifold 'looks the same' i.e. for any 
pair (<I>, of points in V® T there exists an isometry (with 
respect to the Riemannian metric structure) which takes $ into 
>]>. In contrast to spheres the sectional curvature is no longer 
constant, but remains non-negative and bounded from above. 
Since positive curvature contracts distances, the homogeneity 
V® t T suggests certain (with sectional curvature equal to the 
upper bound just mentioned) D-dimensional spheres as simple 
comparison spaces, 



D := di mil ^ T = n t (2T - n t ) 



(17) 



Any code on such a sphere can be transformed into a code on 
the Stiefel manifold without shortening its geodesic distances 
or volumes (compare [11] for details). 

In practice the following has been done. By homogeneity 
in both spaces one 'base' point has been fixed, e.g. the north 
pole N of the sphere and the point ( J ) in the Stiefel manifold. 
All points on the respective manifold can be equivalently 
addressed in the corresponding tangent space (with respect 
to the fixed base point) by utilizing the so-called exponen- 
tial map, which turns radial tangent lines isometrically into 
geodesies emanating from the base point. Representing the 
Stiefel manifold as a U(T) -homogeneous space and similarly 
the sphere as 

S D S 0(D + 1)/0{D) (18) 

(0(D) denoting the group of D x D orthogonal matrices), the 
exponential map exp coincides with the matrix exponential 
exp : A i — ► J2k A k /k\. Restricting cxp to the ball with radius 
7r, exp becomes a diffeomorphism onto the corresponding 
manifolds up to a set of measure zero. Neglecting this sparse 
set, we have now a representation of the sphere and the Stiefel 



manifold in terms of their tangent space at the chosen base 
point. Choosing some orthonormal base {Yi}, resp. {A^} in 
each of tangent spaces induces a canonical mapping from (the 
tangent space of) the sphere to the (tangent space of the) Stiefel 
manifold by 

jvn d9) 

k k 

This mapping is an isometry of tangent spaces and due to 
curvature non-radial distances will be contracted to an extend 
which is determined by the sectional curvature upper bound 
in the Stiefel manifold. Moreover the geodesic distance is 
only an approximation to the 'chordal' distance of the 
receiver. But it can be shown, that both distance measures 
are (locally) equivalent [11]. In summary, spherical codes 
(for example from MWBE sequences in CDMA systems [4]) 
can be transformed into space time codes, whereas spherical 
minimal distances are possibly contracted, but only to an 
extend, which is determined by the (bounded) curvature of 
the Stiefel manifold and the local equivalence of the distance 
measures. 

VI. Lattice spherical space frequency codes and 

DIVERSITY ROTATION 

The construction of high rate codes in systems with large 
block lengths or large number of subcarriers suffers from the 
increasing (brute force) decoding complexity, which becomes 
intractable soon. To fight this, additional structures in the 
code design need to be exploited. One way of embedding 
(algebraic) structure is provided by the design algorithm for 
wrapped spherical codes [5]: This algorithm uses a clever 
and computationally simple one-to-one, distance expanding 
mapping from M. D to the sphere S D C M. D+1 to map lattices 
on (each hemisphere. Having this map and its inverse at 
hand we can describe the overall encoding and decoding 
procedure by Algorithm 1. Note that the lattice encoding 



Algorithm 1 Encoding (1 — ► 8) / Decoding (8 — > 4) 
Require: K no. subcarriers 

n t no. transmit antennas 

L no. channel taps 

ds target (spherical) minimal distance 
Ensure: C space frequency code 

l: D <- nt(2T - nt), T^K/L 

2: Choose some lattice A C ~R D 

3: Scale lattice generator to meet design distance ds 

4: Apply wrapping algorithm in [5] to obtain a spherical code 

(Cs,d s ) 

5: Project Cs to the tangent space TnS d at the north pole 
(orthogonal projection 'upwards', then rescaling to length 
equal to the arclength between the considered point of Cs 
and the north pole N) 

6: Map tangential code points to T^i^V® t>T = T N S D via 

(19); this yields a tangential space time code C t ' ang 
7: Apply the matrix log function to get C' = logC' C 

V n t ,T 

8: Apply (11) with (15), (14) on C', which defines C 



scheme permit large families of code constructions based on 
(apart from the standard cubic lattices Z n ) the root lattices 
(and their duals) A n (A*), D n (D*) in various dimensions 
n, the exceptional lattices #6,7,8. the Coxeter-Todd lattice 
if 12, the Barnes-Wall lattice Ai6, and the Leech lattice A24, 
compare [12] for details. The choice of the design distance 
ds trades minimal distance for rate. For decoding just reverse 
the encoding algorithm until step 4 and apply some lattice 
decoding afterwards. Note that this procedure is suboptimal 
in general due to some points possibly getting lost in step 4 
of the encoding algorithm (compare [13]). But in the high 
rate regime it is expected, that this possible performance 
loss becomes negligible. Furthermore some efforts have been 
made to compensate the 'buffer regions' ([13]) in part: These 
regions originally placed as separating annuli on the sphere 
to guarantee the designed minimal distance ds have been 
reinspected whether formerly rejected 'buffered' code points 
lie at a distance at least ds from all other code points on the 
sphere. 

A common problem of the proposed lattice based con- 
struction is the inherent rigidity of lattices. Although there 
is some control over the minimal distances, there is no direct 
control on the diversity order, i.e. the diversity product p F . 
To overcome the possibility of diversity loss several rotated 
versions of the lattices have been produced as a preprocessing 
step of the encoding algorithm. Of course this procedure 
does not affect the decoding performance. As the axis of 
rotation the 'diagonal' e defined as the sum of all basis vectors 
of the respective lattice has been taken. Denote by W e the 
orthogonal matrix with e\ = W e e whereas ei = (1,0,..., 0)*. 
Set = (0,...,0)* G R D_1 , Ri = (0 °l) and define 
R^- = cxp(aX), a £ R, with X being the antisymmetric 
(D — 1) x (D — 1) matrix with ones on its upper triangular 
part. Then i?i is a rotation of a degrees about the axis ei and 
R(a) = WlR\W e has e as its axis of rotation. 

VII. Simulations and discussion 

For the simulations the channel was modeled for n t = 2 
transmit antennas and L = 2 taps and for simplicity at the 
receiver maximum likelihood decoding has been performed. 
We constructed space time codes with block length T = 4, 
corresponding to lattice dimension D = 12, compare (17). In 
order to obtain full (frequency) diversity for a scenario with 
L = 2 channel taps these codes are designed for systems 
with K = LT = 8 subcarriers. As a reference, an Alamouti 
based space time/frequency code of comparable rate has been 
simulated. Since the block length of the (space time) Alamouti 
scheme is 2 this code family is well adapted to MIMO- 
OFDM scenarios with K = 2L = 4 subcarriers. Other 
choices of K (smaller) would destroy diversity or (larger) add 
redundancy, which lowers the rate and would not give fair 
comparisons. Note that the restriction on K for the Alamouti 
scheme underlines the need for new coding schemes, which 
are capable of meeting more general system requirements 
(K,L,n t ). The spherical space frequency codes proposed in 
this work achieve this by appropriately choosing the spherical 
dimension D. 



Application of Algorithm 1 produced (by the prescription 
of the spherical minimal distance) two sets of spherical space 
frequency codes. The first set consists of very high rate codes 
based on the lattices K12, -D12, -K12 again, and Ayi with 
(space frequency) rates R — 1.95, 1.71, 1.65, 1.64. They will 
be compared to the 8-PSK modulated Alamouti scheme of 
(space frequency) rate R = 1.5. The second set consists of 
medium rate codes based on the lattices if 12, Z12, and A12 
with rates R = 1.32,1.15,1.04 which are compared to the 
QPSK Alamouti scheme of (space frequency) rate R = 1. All 
these lattice have been rotated as described in the previous 
section by an angle a = "/a for a better diversity gain. Figure 
1 shows the simulation results. 
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Fig. 1. Symbol error rate for spherical space frequency codes compared to 
corresponding Alamouti (8-PSK,QPSK) schemes 

The simulations demonstrate that in the high rate regime the 
spherical space frequency codes outperform the corresponding 
Alamouti scheme. Moreover the Alamouti code has the lowest 
(space frequency) rate in this comparison. 

For medium rates the picture changes. Here the Alamouti 
scheme (although the code with the lowest rate in the sim- 
ulation) still outperforms the new spherical space frequency 
codes. Nevertheless, due to the lattice rotation all codes seem 
to have full diversity, even the qubic Z12 lattice performs not 
too bad. 

VIII. Conclusions 

In this work a class of variable rate, fully diverse spherical 
space frequency codes based on the scheme [1] has been 
proposed, which applies to arbitrary channel parameters, i.e. 
no. of subcarriers K, no. of channel taps L, no. of transmit 
antennas n t . The possible rates are determined by the choice 
of the underlying lattice and its designed minimal distance 
(obtained by an appropriate scaling of the lattice). This makes 
the constructions rather flexible to meet a wide range of 
applications. In particular very high rate codes (up to R = 1.95 
for an 8 x 2 space frequency code in the examples) can be 
constructed in this way. The diversity order can be increased 
by some lattice rotation and in the examples full diversity 
has been achieved. The multipath diversity controlled by the 



matrix A (15) has been chosen with respect to an advantageous 
effect on the power distribution among the subcarriers. The 
simulations demonstrate satisfactory performance results (at 
least in the high rate regime), and the lattice based code design 
offers the possibility of moderate decoding complexity. 

However improvements of the code design and its com- 
plexity seem to be possible by careful analysis of the encoding 
algorithm. The connection of space time/frequency codes with 
spherical codes/WBE sequences in CDMA systems comes as a 
surprise, made possible by a differential geometric analysis for 
a (first order) sphere packing approach to space time/frequency 
coding. This offers a new direction in future research efforts. 
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